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Abstract 

We study the centrality dependence to be expected if only charmonium production in 
the corona survives in high energy nuclear collisions, with full suppression in the hot, 
deconfined core. To eliminate cold nuclear matter effects as far as possible, we consider 



g_) ' the ratio of charmonium to open charm production. The centrality dependence of this 

ratio is found to follow a universal geometric form, applicable to both RHIC and LHC in 
collisions at central and forward rapidities. 

> 

' The behavior of charmonium production in high energy nuclear collisions was proposed 

Ti^ij- ■ quite some time ago as a probe of color deconfinement P]. The basic idea was that 

color screening in the hot primary medium prevents the binding of cc pairs produced in 
the early stages of nucleon-nucleon collisions. As a result, the ratio of hidden to open 
charm production should vanish when the energy densities of the medium produced in the 
collision are sufficiently high. Since the observed J/ip production rate is partially due to 
feed-down from higher mass excited states, the suppression should occur in a sequential 
fashion, first for the decay products from ip' and Xc production, then for directly produced 



In subsequent work [5H7J, it was argued that, at the hadronisation point, statistical com- 
bination of the charm quarks present in the deconfined medium could lead to regeneration 
of charmonia. Even if the "primary" charmonium states produced in individual nucleon- 
nucleon collisions are dissociated through color screening, the abundant charm production 
in high energy collisions could provide enough charm quarks to result in "secondary" char- 
monium formation through the binding of combinatorial pairs of c and c quarks produced 
in different initial nucleon-nucleon collisions. Since the charm production rate grows with 
collision energy and centrality, this would imply a corresponding increase in the ratio of 
hidden to open charm. In Fig. [H we schematically compare the statistical recombination 
prediction to that of the sequential color screening scenario. The ratio of hidden to open 
charm in nuclear collisions is normalized to the pp value, scaled by the number of binary 
nucleon-nucleon collisions. The increase of this ratio above unity in central collisions for 
high charm production rates is a crucial prediction of the statistical recombination model. 
On the other hand, the sequential suppression scenario predicts that the ratio should de- 
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crease to values well below the central SPS results at high energy density when direct 
J/ip{lS) suppression sets in. 
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Figure 1: The schematic ratio of J/ip to open charm production in nuclear collisions, 
normalized to the scaled pp value, as a function of the energy density. 

The first RHIC results on charmonium production, after a decade of SPS data, showed 
a rather striking feature. While it seemed clear that the energy density of the produced 
medium should be much higher at RHIC than at the SPS (typical values are 5-6 GeV/fm^ 
at RHIC relative to 2 - 3 GeV/fm^ at the SPS), the magnitude of J/ip suppression in cen- 
tral RHIC collisions was compatible with the SPS data within the measured uncertainties, 
as was the overall centrality dependence. Both the sequential suppression and statistical 
recombination scenarios could explain this behavior with different sets of assumptions: 

• In the sequential suppression approach, full suppression of the excited ip'and Xc states 
is expected, with complete survival of the directly-produced J/ip{lS) states, after 
cold nuclear matter effects are accounted for; 

• The statistical recombination approach predicts a monotonically decreasing survival 
probability for all cc states, coupled with secondary regeneration increasing with 
energy density. 

Predictions for the new data from the LHC in the two scenarios, at still higher energy 
density (at least 8-10 GeV/fm^), lead to relatively distinct outcomes: vanishing mea- 
sured J/ip yield with sequential suppression and considerably enhanced J/ip production 
with secondary regeneration. In both cases, the predictions are cleanest for the ratio of 
charmonium to open charm since, in this case, initial-state effects (shadowing, parton 
energy loss) largely cancel [8]. Because such data are not yet available, the centrality 
dependence of the cold nuclear matter effects introduces some uncertainty. However, the 
presently available data do not rule out the possibility that once all corrections {J/ip feed- 
down from 5-decay, comparable rapidity and transverse momentum ranges) are made, 
the J/ip production rates at the LHC, relative to open charm production, will be quite 
comparable to the corresponding ratios in central RHIC and SPS collisions. 

The aim of this note is to argue that, if indeed the fraction of charmonia survival in nuclear 
collisions is essentially independent of the collision energy, then the survival probability 
at high energy density appears to be a geometric effect due to charmonium production 
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in the corona of the colhsion where the density is too low for dissociation of primary 
charmonium production. 

The so-called corona effect in nuclear collisions has been discussed in quite some detail for 
various observables P4T6]. Its role in heavy quark production has also been considered 
|15] . The basic feature is that, even in a completely central collision, there is a rim region 
in the transverse plane in which the density of wounded nucleons is low, or in which each 
nucleon undergoes very few collisions, perhaps only a single one. In this region, there is 
no produced medium so that, in the corona, all observed effects should be similar to those 
in pp or pA collisions. While in central collisions of heavy nuclei, the ratio of the rim 
area to the total overlap area can be quite small (depending on how the corona region 
is defined), the fraction of the collision region attributable to the corona increases with 
increasing peripherality until, in the most peripheral collisions, there are effectively just 
single nucleon- nucleon interactions (see Fig. [2]). In the aforementioned studies [9HT6], the 
aim was typically to remove or correct for the corona effect. Here we take the opposite 
approach: we assume that the core completely suppresses all charmonium production so 
that only the corona structure and the resulting centrality distribution are relevant. 




Figure 2: Corona contributions for decreasing centrality, as seen in the transverse profiles 
of the colliding nuclei. 



Before proceeding, we briefly recall the definitions of the relevant geometric quantities. 
We use the formulation given in a recent survey ^17j, to which we refer for further details. 
We consider a collision of two nuclei A and B. For a nucleus of mass number A, the 
nucleon density projected onto a plane orthogonal to the collision axis is given by 



Ta{^) = j dzAPA{s,ZA) 



(1) 



Here pA is the probability, normalized to unity, to find a nucleon from nucleus A at the 
location {s,za), where za denotes the longitudinal position and s is the position in the 
transverse plane. For an impact parameter b, the density of the number of collisions in 
the transverse plane is given by 



n,{s, b) = ABrA{s)rB{s - h)a. 
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where is the nucleon-nucleon inelastic cross section. The integral over s then gives 
the total number of coUsions at impact parameter b. Similarly, the wounded nucleon 
density is given by 



n^(s,b) = ATa{s) |1 

+TO(s-b)n 



l-rB(s-b)a^ 
1 - rAis)ai 
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The integral of riyj over s, A^part(b) = / dPsn^is.h), gives the total number of wounded 
nucleons (or participants) for a given centrality b. 

The ultimate magnitude of the ratio of contributions from the corona region to those 
from the overall interaction region depends both on the variables used to distinguish 
between the rim and the central core regions (density of wounded nucleons, collision 
density, energy density) and on the chosen thickness used, i.e., on the definition of the 
corona. Since the radius of a nucleus is generally defined as that distance from the center 
at which the matter density distribution, p^, has decreased by a factor of two, we use 
a similar definition to specify the edge of the corona in nuclear collisions. The resulting 
transverse profiles of the density of wounded nucleons and of the collision density, both 
for central Au+Au collisions, are shown in Fig. |3] for ctj^^ = 42 mb. For this system, 
= nu,{R = 0)/2 = 2.16 fm'^ at R = 5.4 fm while = ndR = 0)/2 = 9.6 fm'^ at 
R = 4.4 fm. 
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Figure 3: Transverse profiles in central Au-Au collisions, showing (a) the density of 
wounded nucleons, n^, and (b) the density of collisions, Uc, in the transverse plane as 
a function of the radius R, calculated with = 42 mb. In both cases, the corona is 
defined as that region in which the density [n^ or Hc respectively) has decreased by a 
factor of two or more. 

A general feature emerges from the geometry of the collision: for a considerable range 
of impact parameters from central through mid-peripheral collisions, the fraction of the 
transerve collision area in the corona relative to the total transverse overlap is small and 
relatively constant. Only in very peripheral collisions, roughly corresponding to 6 > 1.2/2, 
does the fraction increase quite suddenly to unity. To show this more quantitatively, we 
present the ratio of the contribution to the overlap area from the corona relative to the 
total function of collision centrality in Au+Au collisions in Fig. |H Here Fig. Hl^a) 

gives the ratio, R^, of the number of wounded nucleons in the corona to the total number 
of wounded nucleons while Fig. Ht^b) gives the ratio, Rc, of the number of collisions in the 
corona relative to the total number of collisions. In both cases, the centrality dependence 
is shown as function of A^part/A'max) where A^part is the number of wounded nucleons and 

A^max = A^part(& = 0). 
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Figure 4: (a) Ratio of wounded nucleons in the corona region to the total number of 
wounded nucleons. Here the corona region is defined as the region where the density of 
wounded nucleons, n^, has decreased by a factor of two or more relative to the value at 
b = 0, < 2.16 fm~^. (b) Same as (a) for collision density with Uc < 9.6 fm~^. Both 
calculations are for Au+Au collisions with = 42 mb. 



The results shown in Fig. H] employ = 42 mb, appropriate for ^^snn = 200 GeV col- 
lisions at RHIC. For AA collisions at the LHC, a larger cross section is more appropriate. 
We choose = 81 mb. Increasing by nearly a factor of two has only a relatively 
small effect on the behavior of and Rc, as shown in Fig. [51 The value of riuj is indepen- 
dent of cr^ei y see Eq. (jl]). However, the value of R at which nu;(i?)/n^(0) < 0.5 increases 
marginally, causing i?^ to shift to slightly higher values of A'part/A'^max in Fig. [5](a). The 
origin of this shift is thus not due to nw[R) but to the increase in A^part overall since, at 
"^in^ei^ ~^ -^max 2 A whilc finite values of a(^^ reduce A^max relative to this asymptotic 
value. Note that the value of Uc defining the corona is linearly dependent on a^^, see 
Eq. ([2]), so that for ag^^ = 81 mb, the value of = nc(0)/2 is 18.9 fm'^. The value of 
R when Uc = nc{0)/2 is independent of cr^^ . The same shift in A'max described above 
for R^ manifests itself as a steepening of Rc with Apart /Amax for the larger cross section. 
Note that the value of Rc for the two A^A^ cross section values coincide for Apart = A^max. 

The definition of the corona is clearly model dependent and can therefore be given in 
terms of either or Uc- However, charmonium production is a hard process and, as 
such, its rate will be governed by the number of binary collisions. Excluding all other 
possible effects. Fig. Hl^b) can therefore be directly compared to the centrality profiles 
of production data. If we define the corona in terms of wounded nucleons, we have to 
calculate the number of collisions in the corona and core regions defined by in order to 
compare to the data. The resulting ratio, Rcw, of the number of collisions in the corona 
to the total number of collisions is shown in Fig. [6] for both values of (t^^ . 

We note that changing the collision system from Au+Au to Pb+Pb also has only a small 
effect on R^, Rc and Rcw The central values of and Uc will change somewhat. In 
addition, the values of R at which n^/n^(0) < 0.5 and nc/nc{0) < 0.5 are slightly higher 
for the larger A. 
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Figure 6: Ratio of collisions in the corona relative to the total number of collisions for 
Au+Au collisions, using to define the corona. The result is shown for = 42 and 
81 mb. 

Before any comparison to data, we have to address the role of cold nuclear matter effects. 
In the usual treatments, where the charmonium production in nuclear collisions is com- 
pared to that in pp or pA interactions to determine "anomalous" behavior, the role of 
cold nuclear matter effects on charmonium production in both the initial and final state 
has to be taken into account. If one considers instead, as suggested in Ref. [8], the ratio 
of charmonium to open charm, initial-state effects should largely cancel. In this context, 
a first test would be to compare production at central and forward rapidities. At RHIC, 
the AA/pp ratios show differences for the two rapidities. If the hidden to open charm 
ratios are measured as a function of centrality, would the results in the different rapid- 
ity regions coincide? Final-state effects on the charmonium resonances will certainly be 
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reduced in the corona. However, it is not clear if they can be neglected entirely. In partic- 
ular, at the SPS, there are indications that nuclear absorption of the produced charmonia 
in the nuclear medium is not negligible. Analysis of J /if} production at central rapidity 
over the range of fixed-target energies shows that the effective absorption cross section 
decreases with ^/snn This result is not unexpected since, as ^/snn increases, the 
traversal time of the nascent charmonium state through the nuclear medium will become 
very short, diminishing the role of absorption. While absorption could still play a role in 
the comparison between central and forward rapidities, the hidden to open charm rates 
should clarify the situation considerably. 

If we assume that at sufficiently high a/saT/v, the central core of the transverse collision 
region is hot enough to dissolve all charmonium bound states, so that only those produced 
in the corona region survive, then we obtain two immediate predictions for the ratios 
of charmonium to open charm production after the 5-decay contributions have been 
removed. These predictions are the main result of this paper. 

• The J/ ip survival probabilities in central AA collisions for fixed A should be essen- 
tially the same at the SPS, RHIC and the LHC 

• The centrality profiles for J/ip production at RHIC and the LHC should overlap. 
Some differences may be possible for SPS energy if the dissociation in the core is 
incomplete for more peripheral interactions. 

Here we have defined survival as the ratio of ratios: we divide the charmonium to open 
charm production ratio in AA collisions as function of centrality by the corresponding 
single-valued production ratio in pp collisions. No scaling is needed, but if the pp ratio 
is not independent of rapidity, the AA and pp ratios should be compared at the same 
rapidity. 

We can make an additional prediction for the expected transverse momentum distribution. 
Nuclear effects due to initial-state scattering and final-state fiow generally broaden the 
Pt distributions relative to pp interactions. If only corona production remains in high 
energy AA collisions, the px distribution should converge to that in pp interactions. Such 
convergence would discriminate against recombination models of J/ijj production which 
would predict narrower pj- distributions [12]. However, charm fiow has been observed in 
nuclear collisions |20] which, in the case of statistical recombination, should be refiected 
as Pt broadening of the J/ ip produced by recombination. 
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